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Abstract. 

In the search for spin-1/2 kagome antiferromagnets, the mineral volborthite has 
recently been the subject of experimental studies Q] [2| [31 Q| . It has been suggested 
that the magnetic properties of this material are described by a spin 1/2 Heisenberg 
model on the kagome lattice with spatially anisotropic exchange couplings. We report 
on investigations of the Sp(AT) symmetric generalisation of this model in the large AT 
limit. We obtain a detailed description of the dependence of possible ground states 
on the anisotropy and on the spin length S. A fairly rich phase diagram with a 
ferrimagnetic phase, incommensurate phases with and without long range order and a 
decoupled chain phase emerges. 
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Submitted to: J. Phys.: Condens. Matter 

1. Introduction 

The magnetic lattice of the natural antiferromagnet volborthite consists of Cu 2+ ions 
which occupy well separated planar kagome-like nets. A monoclinic distortion leads 
to a slight difference between the exchange couplings along one lattice direction (J) 
and the two other directions (J') within the the planes. In this study, we consider the 
appropriate model, the Heisenberg antiferromagnet on the kagome lattice with spatial 
anisotropy (see figure HJ), 

n, t c c -+i7 J J if i and j ^ c-sublattice 

H = J« S, with J l3 = = . (1) 

n.n. L 

We have applied the spin wave approximation to this model, have done 
quantum perturbation theory around trial groundstates, and have studied the Sp (AO- 
generalisation in the limit of large M where the mean-field approximation becomes exact 

Here, we report the results of the Sp(A/")-symmetric model. We consider this model 
in the full range of the anisotropy < J < oo and pay particular attention to possible 
transitions between different phases of the model which are expected to occur when J 
is varied. 
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Figure 1. Anisotropic kagome model. The coupling J' will be set equal to unity. 8, 
v = 1, 2, 3, are the primitive lattice vectors of the kagome net. 



2. Sp(AT) approach 

In this section, we sketch the mean field treatment of the Sp(A/") generalisation of the 
model (JTJ). The procedure follows closely that of [01 Ej. Further details can also be 
found in j^j. The first step is the replacement of the SU(2) invariant product Sj • Sj 
in ((TJ) by the Sp(A/")-symmetric expression jj>(jT a i3bl a bj)(J'' yS bi 1 bjs), where bi a , b\ a are 
bosonic operators and J is the TV-dimensional generalisaton of the antisymmetric unit 
tensor in two dimensions [01 Ej. Here, we have introduced a paired index notation, 
a = (m, a), m = 1, ■ ■ ■ ,7V, a =|, { for the spinor index of the bosonic operators. The 
boson occupation numbers n b = b\ a b ia are subject to the constraint 

K = n b /M (2) 

with k a fixed parameter For the SU(2) case, J\f — 1, one has k = 2S with 
S the spin length. Hence, large values of k, equivalent to large values of S in the 
SU (2) model, correspond to the classical limit and small values of k correspond to the 
quantum regime. The next step is a decoupling of the quartic terms in the Hamiltonian 
by a Hubbard- Stratonovich transformation using fields Qij which are defined on the 
links of the lattice. Additional fields A have to be introduced as Lagrange multipliers to 
enforce the constraints (j2j) on the sites. We assume that the link variables are identical 
on all horizontal links of figure ^ \Q a b\ = Qi, and on all diagonal links of figure ^ 
\Qac\ = \Qbc\ — Q2- Furthermore, we assume that A = A a on the a and b sites in 
figure [TJ and A = A c on the c sites. Now, the bosons can be integrated out exactly. 
Then, in the large M limit, the physical values of the variables Qi, Q 2 , A a , and A c are 
determined by the saddlepoint of the mean field energy 

= J|g 1 | 2 + 2|g 2 | 2 -(2A a + A c )( K + l) + -^^^(k)(l + Mk)| 2 ) (3) 
yviV v iV v M 

with respect to these variables. Here, iV v is the number of downward pointing 
triangles of the lattice. The spinon frequencies ^(k), /j, = 1, 2, 3, are obtained as 
the three positive solutions of detD(o>) = 0, where D(a>) is the six dimensional matrix 
related to the bosonic part of the Sp(jV) Hamiltonian after the Hubbard-Stratonovich 
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transformation has been performed: 
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with <3 M (k) = iQ^ sin(<5 At k/2), /i = 1, 2, 3, and Q% = Qi- In the above formula, 
the matrix elements Q M (k) incorporate information on the lattice symmetries via the 
primitive lattice vectors <5 M of the kagome net, see figure and the wave vector k. The 
quantity ^(k)! 2 represents the density of Bose condensate that will develop at those 
points in the Brillouin zone where one of the spinon frequencies vanishes. By calculating 
the Sp(7V) analog of the spin-spin correlation function, we can show that the long 
distance behaviour of this function is directly related to the existence of a condensate 
[Hj: a finite condensate density, |x M (k)| 2 > implies long-range order (LRO), while for 
^(k)! 2 = there is only short-range order (SRO). 

To find the saddle point of Emf, ©, in the space of the variables Q\, Q2, A a , 
A c and |x M (k)| 2 at a given pair of parameters J, k, we use an iterative numerical 
routine. Remarkably, and in contrast to previous applications of the mean field Sp(AT) 
approach [HI IE], it turns out to be indispensable to allow for two independent Lagrange 
multipliers A a and A c . In view of the anisotropy of our model this is plausible. 



3. Results and discussion 



The results of the large-A/" approach are summarised in the zero temperature phase 
diagram of the anisotropic kagome antiferromagnet. It displays the phases that we were 
able to discern in the J — 1/k plane. We first note that as expected LRO which is always 
present for large k disappears from the incommensurate (IC) phase for sufficiently small 
values of n. Recall that for large k we are in the classical regime of our model. The 
phase boundary that separates the region with SRO from the region with LRO was 
found by determining for each value of J that value of k where the gap in the lowest 
branch of the one spinon spectrum ^(k) first vanishes as k is increased. At larger 
values of k, condensate, and hence LRO will appear. Since our model is maximally 
frustrated for J — 1, one expects a maximal suppression of LRO by quantum effects at 
this point. As is seen in figure El this is confirmed in the large- M approach. For J = 0, 
the exact quantum ground state of our model ([TJ is ferrimagnetic (FM) according to 
the Lieb-Mattis theorem pj. The classical ferrimagnet has all spins on the horizontal 
lattice lines aligned ferromagnetically while the middle spins are oriented antiparallel 
to these. In the classical limit, k — » 00, of our Sp(A/")-model, the expectation value 
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Figure 2. Phase diagram; FM: ferrimagnetic phase, IC: incommensurate phase, DC: 
decoupled-chain phase; empty circles, full circles, and asterisks: numerical results; full 
lines: analytical results. Asterisks mark the boundary between phases with long-range 
order (LRO) and short-range order (SRO). 



Qi = {J ma m> a ' b\ ma bX m i ■, which measures the singlet weight on the horizontal bonds, 
should vanish in this state for sufficiently small J, J < Jfm( k )- As J is increased 
beyond Jfm, Qi increases in the manner of an order parameter at a second order phase 
transition (see panel (a) of figure EJ) • At the same time, the parameter Q 2 begins to 
decrease, and it drops to zero at J = Jdc( k ) ( see panel (b) of figure EJ). (Note that the 
J-scales are different in panel (a) and panel (b)). Thus, the large- M approach predicts 
the existence of a decoupled-chain phase in the region above the phase boundary J D( > 
Q2 decreases to zero continuously so that the phase transition at Jdc appears to be of 
second order again. Both, LRO and SRO phases may be characterised by an ordering 
wave vector q or d = 2k min , where k min is that wave vector at which the lowest branch of 
the one-spinon excitation spectrum ^(k) has its minimum. The static spin structure 
factor S(q) develops a peak at q or d- We display q ord in panel (c) of figure El For 
J < Jfm(^) we find q* rd = —2tt which is equivalent to q ord = 0. For J > Jfm(^), 
q x OTd increases with J, and it reaches the value — 47r/3 at the isotropic point J = 1. 
There it is independent of the value of k in agreement with Ref. [Zj. For k > 1/3, the 
behaviour of q x ord as a function of J is as expected: it increases monotonously towards 
\lord\ = vr as J ^ 00. However, for n < 1/3, the function ql rd (J) develops a minimum 
at J m 1.5 which becomes more pronounced as k decreases, i. e. as the importance of 
quantum effects increases. At present we have no physical explanation for this behaviour 
of q or d- The behaviour of the saddlepoint values of the Lagrange multipliers A a and A c 
as functions J is displayed in panel (d) of figure El for various values of k. As we have 
emphasised above, A a 7^ A c for general anisotropy, J / 1. Only at the isotropic point 
we find A a (ft) = A c (k) in accordance with the expectation. Again, we have no intuitive 
physical explanation for the behaviour of these parameters as functions of J and k. 
Quite unexpectedly A c drops to zero at the boundary Jdc of the DC phase. 

In the vicinity of the FM-IC and the IC-DC phase boundary, where Qi respectively 
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Figure 3. Saddlepoint values of the parameters Q±, Q2, of the ordering wave number 
Qord (lord ~ 0)> an< ^ °f ^° an< ^ Full circles and horizontal lines in panels (a), (b), 
and (c): results of analytical determination of the saddlepoint. 



Q 2 , are small one can expand E MF , ©, with respect to these variables and solve the 
extremum problem analytically. The expansion reads 

EuFiQa) = e 0a + r a \Q a \ 2 + g a |Q a | 4 + 0{\{Q a \ 6 ) , (5) 

where a = 1(2) for the FM-IC(IC-DC) phase boundary. After the stationarity conditions 
Emf(Qo) with respect to the variables Qp, [3 = 2(1), A a , A c have been solved the 
coefficients r 0a and g a are known as functions of J and k. Our numerical results suggest 
that g a > 0. Similarly as in a Landau- Ginzburg expansion, the FM-IC and the IC-DC 
phase boundaries are then obtained as the solutions of r aQ = 0. These solutions are 
presented as solid lines in figure 121 Obviously, the numerical data agree reasonably well 
with these analytic solutions. Small discrepancies arise in the numerical data because 
the energy surface E MF (Qi, Q 2 , X a , A c ) is very shallow in the vicinity of the saddle point. 
On these boundaries and inside the FM and the DC phases the values of Q\, Q2 and 
of the ordering wave vector remain constant as is shown in the corresponding panels 
of figure H3 Using a similar expansion technique for the Sp (Af) -analog of the spin-spin 
correlation function, we have also been able to obtain the asymptotic behaviour of this 
function on the phase boundaries analytically. We find that there is indeed LRO along 
the entire phase boundaries up to arbitarily small values of k. However, while along the 
FM-IC phase boundary there is LRO between any pair of spins S^j, S u j, fj,, v — a, b, c, 
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along the IC-DC phase boundary, there is LRO only between the spins on the c sites 
S c j, see figure [TJ 

To conclude we wish to compare the results presented here for the anisotropic 
kagome antiferromagnet (AF), (JTJ), with those obtained by Chung et al. [Hj for the Sp(AT) 
symmetric generalisation of the anisotropic triangular AF. With increasing values of J, 
which controls the anisotropy in both models, we find the same succession of phases as 
in ref. [H]: a collinear phase, an incommensurate phase and a decoupled chain phase. 
However, while LRO disappears for sufficiently small values of k = 2S for all values of J 
in the triangular AF, we find LRO in the entire FM phase and along the FM-IC and the 
IC-DC phase boundaries up to arbitrarily small values of k. For both models, the mean- 
field Sp(JV) approach predicts a DC phase which, however, is certainly not described 
faithfully by this approach. Qualitative consideration of the role of fluctuations in this 
phase lead the authors of ref. [Hj to the conclusion that there is spin Peierls order in 
the DC regime of the anisotropic triangular AF. Such a picture is still elusive for the 
anisotropic kagome AF. However, by applying perturbation theory to the trimerised 
version of the model (0), i. e. to a model for which the couplings on e. g. the upward 
pointing triangles are uniformly weaker than those on the downward pointing triangles 
[TU] , we find that for J ^> 1 the dynamics of the c spins of our model, see figure ^ is 
effectively described by the Hamiltonian of the anisotropic triangular AF with strong 
coupling in one of the three lattice directions. Therefore, according to the arguments 
of ref. [S] there should be spin Peierls order between these spins. Since, on account of 
the trimerisation, the spins on the a and b sites in figure Q are already paired in nearest 
neighbour singlet pairs, one may thus speculate that spin Peierls order also prevails in 
the DC phase of the anisotropic kagome AF. 
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